Spectral properties of a spin-incoherent Luttinger Liquid 
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We present time-dependent density matrix renormalization group (DMRG) results for strongly 
interacting one dimensional fermionic systems at finite temperature. When interactions are strong 
the characteristic spin energy can be greatly suppressed relative to the characteristic charge energy, 
allowing for the possibility of spin-incoherent Luttinger liquid physics when the temperature is high 
compared to the spin energy, but small compared to the charge energy. Using DMRG we compute 
the spectral properties of the t — J model at arbitrary temperatures with respect to both spin and 
charge energies. We study the full crossover from the Luttinger liquid regime to the spin-incoherent 
regime, focusing on small J/t, where the signatures of spin-incoherent behavior are more manifest. 
Our method allows us to access the analytically intractable regime where temperature is of the 
order of the spin energy, T ~ J. Our results should be helpful in the interpretation of experiments 
that may be in the crossover regime, T ~ J, and apply to one-dimensional cold atomic gases where 
finite-temperature effects are appreciable. The technique may also be used to guide the development 
of analytical approximations for the crossover regime. 



PACS numbers: 71.10.Pm,71.10.Fd,71.15.Qe 



I. INTRODUCTION 



A remarkable aspect of one-dimensional interacting 
electron systems (we will use one-dimensional electrons 
as a concrete example throughout the paper, but our 
results immediately generalize to any system with an in- 
ternal "spin" degree of freedom, e.g. cold atomic gases) 
is that a perturbative treatment of the interactions about 
the non-interacting limit, and thus Fermi liquid theory, 
fails. This is a consequence of the pervasive nesting tak- 
ing place at all densities and polarizations, due to the 
simple fact that the Fermi surface reduces to two Fermi 
points. The central result to emerge in one dimension 
is that, for many realistic parameters, the low-energy 
physics is gapless and described by a universal low-energy 
theory called "Luttinger liquid" (LL) theory^^ Accord- 
ing to LL theory, there are no electron-like quasi-particles 
analogous to those found in Fermi liquid theory (which 
describes interacting electrons in three dimensions). In- 
stead, the low energy physics is dominated by bosonic 
collective excitations. LL theory also states that for fi- 
nite interactions there will be a spin-charge separation 
with distinct collective spin and charge excitations that 
each have their own characteristic velocity and Hamilto- 
nian. The spectral properties of the LL are very different 
from a Fermi liquid, but have been computed and are 
known 

A particularly good realization (low disorder) of one- 
dimensional electrons is found in high mobility semicon- 
ductor heterostructures of the type often used to study 
the fractional quantum Hall effect. Related systems were 
used recently to establish the presence of LL physics 
in quantum wires^Si While previous carbon nanotube 
experiments^^ demonstrating a power-law form of the 



tunneling density of states were correctly interpreted as 
an indication of LL behavior, they did not unambigu- 
ously establish its existence because they did not probe 
the full spectral function of the system due to the lo- 
cal tunneling of electrons (which docs not allow momen- 
tum resolution). The key feature of the semiconductor 
heterostructure devices is that parallel wires can be fab- 
ricated and momentum-resolved tunneling experiments 
performed. It is the momentum resolution that allowed 
the dynamical properties of the wires to be measured and 
LL behavior to be unambiguously observedi^ 

However, these experiments also showed a distinct set 
of behaviors when the temperature was estimated to 
be of the order or much larger than the characteris- 
tic spin energy^ In this regime, LL theory is not ex- 
pected to hold, but rather a separate theory describ- 
ing "spin-incoherent" electrons takes overJ^ It turns out 
this "spin-incoherent Luttinger liquid" (SILL) has many 
more universal properties than the LL^ but its conclu- 
sive demonstration in experiment is not yet universally 
agreed upon^ Part of the challenge is that for realis- 
tic parameters many one-dimensional systems fall in the 
crossover regime between LL and SILL, greatly compli- 
cating the interpretation of the experiments! 14 ' 15 This 
crossover regime is not easily or accurately handled by 
existing analytical methods, so a numerical approach is 
required. 

In this work, we describe a technique well suited to 
this challenge and compute several quantities that can 
be directly compared to experiment. The qualitative 
agreement with existing experiments in Ref . [Ill ] provides 
further evidence that the spin-incoherent regime has in- 
deed been reached. At present, there are no other nu- 
merical methods that have been demonstrated to accu- 
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FIG. 1: Momentum resolved spectrum at zero temperature 
of a t — J chain of length L — 64, with N = 48 particles, 
and (a) J = 0.5 (b) J = 0.05 (in units of t), obtained with 
time-dependent DMRG. Negative frequencies correspond to 
the photoemission spectrum obtained by removing a fermion, 
while positive values correspond to inverse-photoemission. 
The spinon bands have a weak dispersion of width ~ J. 
Holon, and corresponding shadow bands, are clearly visible. 
Frequencies are measured in units of t. 



II. SPIN-INCOHERENT HUBBARD CHAIN 

To be concrete, we study the one-dimensional Hubbard 
model: 

L L 

H = -t^2 ( c L c m<7 +h.c.J +[/^V; T n u , (1) 



where c it7 creates an electron of spin a on the i th site 
along a chain of length L. The hopping parameter of the 
Hubbard chain is t, the onsite interaction energy is U, 
and we take the inter-atomic distance as unity. In the 
limit of large repulsive U, we can cquivalcntly consider 
the t — J model, defined as 
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h.c.) 
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(2) 
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where the constraint forbidding double-occupancy has 
been imposed. The natural excitations of this model are 
charge and spin collective modes (holons and spinons, 
respectively) with different velocities that depend on the 
ratio U/t, or J/t. In the U — > oo, J — > limit, the 
ground state factorizes into the product of a fermionic 
wave function \<f>), and a spin wave function |x)^ 



|g-s-. 



® lx>- 



(4) 



The first piece, \<f>), describes the charge degrees of free- 
dom, and is simply the ground state of a spinless non- 
interacting tight-binding Hamiltonian. At finite U , the 
spins are governed by a Heisenberg interaction 



H, 



Si+i 



(5) 



rately access the parameter regimes and system sizes we 
study here. A strength of our density matrix renor- 
malization group (DMRG)-based calculation is that it 
free of the artifacts introduced by statistical sampling, 
such as in quantum Monte Carlo. Our calculations start 
deep within the spin-incoherent regime and approach the 
crossover regime "from above" (i.e. from temperatures 
above than the crossover temperature), complementing 
existing analytical methods that attempt to approach the 
crossover regime "from below" using LL theory. Most im- 
portantly, the method accurately captures the crossover 
regime allowing its properties to be revealed. Such in- 
formation can be crucial in the proper interpretation of 
experiments in the strongly interacting regime (which are 
always at finite temperature) and can be used as an aid 
in the development of approximate analytical methods to 
describe this regime. 



where J depends on the charge wave-function and is pro- 
portional to 4t 2 /UiH In the U — > oo, J — > limit, the 
spin states are degenerate and the dispersion is just a 
non-interacting band e(k) = — 2icos(fc), but any finite 
interaction will lift this degeneracy and give the spin de- 
gree of freedom some dispersion. The factorized wave- 
function approach has been used in a number of key ear- 
lier studie o 16 ' 17 i 20 ' 21 in the strongly interacting regime, 
but we do not make such an approximation here. 

In Fig[T]we show the momentum resolved spectrum of 
a chain with L = 64 sites and N = 48 particles, for values 
of J = 0.5 and J = 0.05 (all quantities are in units of 
the hopping), obtained with time-dependent DMRG at 
zero temperaturej 18 i 19 The spectrum is clearly gapless, 
displaying a weakly dispersive spinon band of width ~ 
J, and broad holon bands of width ~ it. Our results 
agree in the U —> oo limit with the dispersion calculated 
m Refs.^H using the factorized wave function, Eq. 
((31); and also with the exact diagonalization results for 
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FIG. 2: Specific heat of a t — J chain of length L = 32 and 
N = 24 fermions, calculated with time-dependent DMRG. 
Temperature is in units of the hopping t. 



the t — J model in Ref. [20|, at the same density and 
similar value of parameters. For small values of J we 
see an almost non-dispersive spinon band. In this case, 
a spin-incoherent behavior would be observed at a finite 
temperature larger than the characteristic spin energy 
scale, but much smaller than the Fermi energy J<T« 
Ep ~ t. When these conditions are realized, the spins 
are totally incoherent, effectively at infinite temperature, 
while the charge sector remains very close to the ground 
state. 
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FIG. 3: (a) Spin structure factor, and (b) momentum distri- 
bution of a t — J chain of length L = 32 and N = 24 fermions, 
with J = 0.05, for different values of the temperature. The 
thick full lines in (a) and (b) correspond to T = 0. Arrows in- 
dicate increasing (3 (decreasing temperature) in units of 1/t, 
in steps of four, (c) shows the behavior of the "Fermi mo- 
mentum" k* F as a function of temperature, obtained as the 
inflection point in the momentum distributions shown (b). 



III. METHOD 

The key idea behind our calculation is thermo field 
dynamics . 22 i 23 ' 24 i 25 ' 26 i 27 ' 28 This construction allows one 
to represent a mixed state of a quantum system as a 
pure state in an enlarged Hilbcrt space. Consider the en- 
ergy eigenstates of the system in question {n}, described 
by a Hamiltonian H , and introduce an auxiliary set of 
fictitious states {n} in one-to-one correspondence with 
{n} . We can then define the unnormalized pure quan- 
tum state, 

|^G9)) = e-^XO)) = £ e-^' 2 \nn) (6) 

n 

where n is a copy of n in the auxiliary Hilbert space, (3 = 
1/T is the inverse temperature, and \ip(0)} = ^2 n \nn) 
is our thermal vacuum. Then the exact thermodynamic 
average of an operator O (acting only on the real states) , 
is given by 

(d) = z(p)- i {m\d\m), (7) 

Where the partition function is the norm of the thermal 
state Z{(3) = {^{(3)\tp{(3)) . We can clearly see how the 
calculation of a thermodynamic average reduces to calcu- 
lating a conventional expectation value of an operator in 
a pure quantum state, at the price of working in a larger 
Hilbcrt space. 

At j3 = 0, the state | -0(0)} is the maximally entangled 
state between the real system and the fictitious system. 



We can see that this is independent of the representation, 
and we can choose any arbitrary basis. In particular, is 
natural to work in an occupation number representation 
where the state of each site i takes on a definite value rt;. 
One finds 

iv>(o)>=nEi n ^>=ni j °*>' ( g ) 

i rii i 

defining the maximally entangled state \Ioi) of site i 
with its "ancilla", the local degree of freedom in the 
auxiliary system. At this point it becomes convenient 
to perform a time-reversal transformation on the an- 
cillas. Therefore, for the case that concerns us, where 
double-occupancy is forbidden, this state can be written 
as \I 0i ) = | t,I> - | -l,T) + |0,6). This simple step al- 
lows us to work in a basis where the total spin projection 
S* ot of the chain-ancilla system is effectively zeroi 2 ^ We 
emphasize that both spin and charge degrees of freedom 
appear in \Ioi), and are therefore treated on equal footing 
as regards finite-temperature effects. (Subject, of course, 
to the no double occupancy constraint.) 

The state of the system at an arbitrary temperature 
(3 is obtained by evolving the maximally mixed state in 
imaginary time, Eq. ([5]) with (3 = 0, using the Hamilto- 
nian acting on the real degrees of freedom. The ancillas 
do not have any interactions controlling their dynamics. 
They evolve only by their entanglement with the phys- 
ical spins, effectively acting as a thermal bath. This is 
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the basis of the finite-temperature DMRG method^ No- 
tice that at zero temperature, the site and the ancilla are 
totally disentangled, while at finite temperature there is 
always a finite degree of entanglement that only depends 
on the dynamics of the system. 

An important consequence of the previous description 
is that it would correspond to working in the grand 
canonical ensemble: Even though the spin and charge 
quantum numbers are conserved for the enlarged system, 
this is not the case if we restrict ourselves to the physi- 
cal chain. In order to work in the canonical ensemble we 
need to start from a thermal vacuum where the physical 
states \n) and their copies \h) have each a fixed number 
of particles. To achieve that, we are going to construct 
a state that is a sum of all possible states of charge and 
spin, with the constraint that the total number of par- 
ticles on the chain has to be equal to N , and that the 
charge state of the ancillas is an exact copy of the charge 
state of the physical chain. We achieve this by calcu- 
lating the ground state of a very peculiar Hamiltonian, 
using conventional DMRG: 

H = -^(AjA J+ h.c). (9) 

The operator A T (A) creates (annihilates) a singlet be- 
tween the physical spin and the ancilla, 

A ! = (4A - 4A) m ( 10 ) 

where the "tilde" operators act on the ancillas on site 
i. The ground state of this Hamiltonian is precisely the 
equal superposition of all the configurations of N "phys- 
ical site-ancilla" singlets on L sites. This state can be 
represented very efficiently in terms of a matrix-product 
state, and consequently, by the DMRG method. In prac- 
tice, the number of DMRG states required is of the order 
of the number of particles. We find the use of the "entan- 
gler" Hamiltonian practical and convenient. Note that it 
does not disrupt the SU(2) symmetry of the t— J model. 

IV. GREEN'S FUNCTIONS 

We study the spectral properties of a spin-incoherent 
chain by evaluating the Green's functions at time t at 
finite spin temperature: 

G(x-x o ,t,0) = (i } (p)\e lHt -' t dHx)e~ lHt -' t d(x )\^((3)), 

where the generic operators of interest 0(x),O' ! (x) act on 
the system at site x, and the Hamiltonian H t -j governs 
the physics of the actual physical chain, not including the 
ancillas. 

We use a similar method to the one described in 
Ref. {30l| . The calculation proceeds as follows: First, we 
evolve the maximally entangled state in imaginary time 
to the desired value of (3 measured in units of 1/t, e.g. 



(3 = 2 means T = t/2. Then, an operator 0{xq = L/2) 
is applied in the center of the chain. The resulting state 
is evolved in real-time, and at every step we measure the 
overlap with the state O(x)e~ lHt - jt \ip(0)) . We obtain 
the desired Green's function in frequency and momen- 
tum by Fourier transforming the results in real space and 
time. Both states, \ip(f3)), and O(xo)\ip(f3)), have to be 
evolved in real time. In this work we use a third order 
Suzuki- Trotter decomposition with a typical time-step 
t = 0.1, both for the real-time and imaginary-time parts 
of the simulation, keeping 800 states, enough to main- 
tain the truncation error below 10~ 5 . As customary in 
most DMRG calculations, we used open boundary condi- 
tions, and by doing the Fourier transform we are assum- 
ing that boundary effects can be ignored, as though the 
system were translational invariant. In order to minimize 
the finite-size effects induced by the boundarie o 31 ' 32 ! 33 
we evolve to times t = 15, and Fourier transform to fre- 
quency using a Gaussian window or width a = 6 in the 
time domain, which in turn leads to a mode getting an 
artificial broadening in frequency proportional to 1/er. 
We point out that we have not used the linear predic- 
tion method introduced in Ref. [30], but the bare data 
obtained from the simulation. At zero temperature we 
have found that this works well, reproducing the features 
observed in the Bethe Ansatz solution of the Hubbard 
chain (compare our Fig. 2(c) to Fig. 7 in Ref. [3), namely, 
the singularity at 3/cf (seen at 2ix — 3kp). At finite tem- 
peratures the system develops a finite correlation length, 
which is further enhanced at higher temperatures due to 
the spin-incoherent mechanism - to be discussed below. 
This is reflected in a localization that makes the bound- 
ary effects irrelevant. Working with open boundary con- 
ditions also avoids the degeneracy occurring in systems 
with periodic boundary conditions and size L = An, with 
n being an integer.— The numerical errors can be at- 
tributed to the accumulation of truncation error, and 
the Trotter decomposition. The latter are under control, 
while the truncation error would translate into error bars 
that are much smaller than the broadening in frequency, 
and are therefore ignored for visualization purposes. 

V. RESULTS 

In Figs [2] and [3] we show some characteristic physical 
quantities at finite temperature, such as the specific heat 
Cy, spin structure factor S{k), and momentum distri- 
bution function n(k), for a chain with L = 32 sites and 
N = 24 fermions (3/8-filling). All results correspond to a 
value of J = 0.05. The correlation functions are defined 
as: 

S(k) = ^(S!S^}e ik ^, (12) 
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FIG. 4: Photoemission spectrum of a one-dimensional t — J 
chain with L — 32 sites and N = 24 fermions, and J = 0.05. 
Different panels show different values of temperature T = 
1/13. The crossover to the spin-incoherent regime is achieved 
at f3 ~ 20. 



The specific heat in Fig[2] shows a clear peak at a value 
of the temperature T ~ J, signaling the onset of the 
spin-incoherent regime where the spin degrees of free- 
dom are highly thermalized. At larger temperatures, the 
spin degrees of freedom are saturated, but a broad peak 
associated with the charge degrees of freedom is appar- 
ent. The spin structure factor in Fig[3Ja) shows a peak 
at momentum k = 2fcp that develops precisely at low 
temperatures T < J. The zero temperature result is 
consistent with that computed in Ref.pjj]. The momen- 
tum distribution in Fig|3fb,c) displays the expected Lut- 
tinger liquid profile, with no discontinuities at the Fermi 
point. Below this value of the temperature we also notice 
the onset of a singularity at k = 3kp in the momentum 
distribution* 3 ^ This singularity corresponds to the trans- 
fer of spectral weight to the shadow bands, that origi- 
nate from the scattering with the spin fluctuations that 
diverge at k = 2kp i 16 i 20 i 21 While this behavior had al- 
ready been seen in finite-temperature calculations using 
the factorized wave function ^ in Ref.[36[, our calcula- 
tions do not rely on the factorized wavefunction, or the 
XY approximation in the spin sector (rather than full 
Heisenberg symmetry). Moreover, they can be readily 
generalized to a number of other spin symmetries, in- 
cluding the incorporation of spin-orbit effects. 



We approximated the temperature- dependent Fermi 
momentum k* F by taking the inflection point where n{k) 
changes curvature, and plotted the result in FigEKc). 
The Fermi momentum moves continuously from the 
zero-temperature value kp = irN/2L to 2kp, with the 
crossover region centered around T ~ J, as expected^! 
We actually observe a saturation value below 2kp, but 
this is an artifact of taking the non-rigorous definition 
of kp as the inflection point in n(k). It is important to 
note that n(k) for the Hubbard or t — J models changes 
its form qualitatively for fillings larger than 1/4, but less 
than 1/2 .±2 

The "special" value of 1/4 filling in the lattice mod- 
els (as opposed to the effectively low-energy theories) is 
related to the underlying Hubbard model, and the qual- 
itative change for fillings above it and below it can be 
understood in the following way: The shadow bands 
carry a significant ammount of spectral weight above 
k = kp. When the density is larger than quarter- filling, 
the shadow band covers all fc-space from to it. This 
translates into and increase of weight above kp. (Recall 
that n(k) is the integrated weight for uj < 0i 20 i 21 ' 35 ) 

One further feature of n(k) is particularly striking: The 
values n(kp) and n(2kp) are temperature independent 
within the accuracy of our calculation, with kp = ttN/2L 
and 2kp twice rather than the value obtained by the 
inflection point of n(k). The same behavior was also ob- 
served in the factorized wavefunction approach with XY 
symmetry described in Ref.[36]. Evidently, then, it is 
independent of the spin symmetry, and probably results 
from an effective temperature independence of the charge 
sector (temperature is taken to be precisely zero with re- 
spect to the charge sector in Ref . [3a] ) . Since both calcu- 
lations are effectively in the large U limit, it may be that 
the temperature independence of these two points can be 
attributed to two extreme "charge configurations" -one 
"evenly spaced" and one " maximally paired" (two sitting 
right next to each other). In the density-density corre- 
lations, the former would correspond to Akp oscillations 
and the latter to 2kp oscillations. In the large U limit, it 
must be that these are both extreme "spin-independent" 
configurations in the sense that "evenly spaced" electrons 
or "maximally paired" paired electrons have only mini- 
mal contributions from the spin energy, leading to the 
temperature independence of n(kp) and n(2kp). 

Qualitatively, the shift from kp to 2kp (as measured 
by the inflection point of n(k)) when the spin-incoherent 
regime is obtained can be understood as a shift from 
particles with spin dynamics to particles that are effec- 
tively spinlessi 37 i 38 In the large but finite U limit of the 
Hubbard model, electrons at zero temperature "dimer- 
ize" ever so slightly and in this way maintain a "mem- 
ory" of their non-interacting kp. However, once T ~ J, 
this dimcrization is washed out (because the energy scale 
for dimcrization is set by J) and effectively shifts kp to 
its "spinless" value, 2kp.— 

Fig0] shows the momentum resolved photoemission 
spectrum obtained by taking 0(x) — c^(x) in Eg. pip 
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FIG. 5: Photomoemission results for the same system as 
FigfJ] and /3 — 10. The right panel shows the same data 
with the amplitude represented as lines. Note the qualitative 
agreement with the experimental results of Ref.[llfl. This 
indicates that disorder effects are not needed to explain the 
main features of the data, and spin-incoherent physics is likely 
relevant. 

in the previous treatment. At infinite temperature, it 
resembles a band of non-interacting spinless fermions, 
following a — 2icos(fc) dispersion with a "width" much 
larger than seen in a zero temperature calculation {e.g., 
the result in FigQ]for a larger system size). As the tem- 
perature is lowered, and (3 increased, we see spectral 
weight being transfered from positive to negative ener- 
gies. At the same time, the band appears to broaden 
in the momentum direction, also splitting into seemingly 
discrete weakly dispersive levels. At a value of (3 ~ 10, 
the strongest features of the dispersion describe a number 
of discrete levels on top of a band that follows the spinless 
dispersion, with very small weight above the Fermi level. 
The high weight at zero momentum corresponds to the 
high density of states in the van Hove singularity of the 
spinless band. At temperatures below T = 1/20 = J, the 
dispersion splits into two "echoes" , centered at k = ±kp, 
showing the emergence of the shadow bands and features 
more reminiscent of a LL. We have verified that the gap 
between the horizontal levels are a finite-size effect, and 
the spacing grows as 1/L as we reduce the system size. 
The two-peaked features for (3 < 10 in the horizontal 
dispersion correspond to scattering of charge states with 
the non-dispersive spins present in the spin-incoherent 
regime. 

Many of these features can be qualitatively understood 
within SILL theory. First imagine a system at zero tem- 
perature with J < t. At the fillings we consider, this 
system will behave as a LL because it is gaplcss. The 
spectral function will exhibit cusp-like singularities at 
lj = v a k and uj = v p k where v a is the spin velocity and 
v p is the charge velocity 4£ Since J « t, »„ < »,,. The 
"sharpness" of the cusps are determined by the interac- 
tion parameters of the LL theory in the spin and charge 
sectors.^ If one now adds a small finite temperature (so 
as to remain in the LL regime) the LL correlation func- 
tions obtain a finite correlation length £ oc v a /T. This 
correlation length will smear and broaden the cusp-like 



singularities £ As the temperature is further raised, there 
is a smallest correlation length than can be obtained: the 
interparticle spacing. In the spin-incoherent regime, £ ef- 
fectively saturates at this value and leads to a universal 
broadenin g 38 i 39 ' 40 of ~ In (2) fop /V of the singularity as- 
sociated with the charge mode and a vanishing of the 
singularity associated with the spin mode. This effect is 
evident in Fig f5] for (3 = 1 when one compares to the zero 
temperature result in Fig[TJ For (3 ~ l/J the shadow 
bands are beginning to emerge and the spin degrees of 
freedom are starting to become dynamical leading to a 
complicated spectral form. 

Finally, we focus on the crossover regime at [3 = 10. 
The spectrum is shown with clarity in Figj5] These re- 
sults can be compared to experiments in nanowiresJ^ 
(See e.g., their Fig. 3.) The qualitative agreement indi- 
cates that the experiments were most likely in the regime 
of highly thcrmalizcd spin states. Moreover, our calcula- 
tions conclusively demonstrate that disorder effects need 
not be invoked to explain the data. n i 13 These results can 
also be compared to high-energy angle-resolved photoc- 
mission experiments on quasi one-dimensional SrCu02, 
where the V-shaped dispersion is also observed4i 



VI. SUMMARY AND CONCLUSIONS 

We have presented a numerical study of the spec- 
tral properties of t — J chains at finite temperature, 
using a generalization of time-dependent DMRG tech- 
niques that combines evolution in real and imaginary 
time. The study of finite temperature effects on the spec- 
tral functions of one-dimensional systems using quan- 
tum Monte Carlo technique a 41 i 42 ' 43 has mostly focused 
on the interpretation of photoemission experiments on 
quasi one-dimensional compounds such as SrCu0 2 44 ' 45l4s 
and TTF-TCNQ41 While the Monte Carlo technique is 
free of the sign problem in one-dimension, the calculation 
of spectral properties involves an analytic continuation 
from Matsubara frequencies which is not straightforward 
in the spin-incoherent regime. The application of a max- 
imum entropy method to the results is affected by sta- 
tistical uncertainties, inherent from the stochastic QMC 
approach. On the other hand, our method is naturally 
applied to study the spin-incoherent regime and we have 
demonstrated it is quantitatively accurate by comparison 
with Bethc ansatz results in various limits. 

We have clearly seen that at temperatures of the order 
of T ~ J, the system experiences a crossover from a spin- 
coherent to a spin-incoherent regime, which is clearly 
manifest in the spectra. Our results in finite systems 
show a compelling qualitative agreement with experi- 
ments in nanowiresAi The fact that our systems have 
a finite size works to our advantage since our parameters 
are similar to the experimental conditions, which involve 
wires at low densities with few electrons. 

In summary, we have been able to address an impor- 
tant and analytically inaccessible regime of strongly cor- 
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related one-dimensional systems. The time-dependent 
DMRG method has the power to access the full crossover 
from SILL to LL behavior as a function of temperature 
and therefore is a powerful tool in the interpretation of 
experimental results, and as a guide to analytical approx- 
imations not yet developed. The technique can be readily 
adapted to study a number of related problems, includ- 
ing cold atomic gases which are notoriously plagued by 
finite temperature effects 4£ 
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